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Abstract. Quantum electrodynamics (QED) fixed in the 't Hooft-Veltman gauge is renormal- 
ized to three loops in the MS scheme. The /3-functions and anomalous dimensions are computed 
as functions of the usual QED coupling and the additional coupling, ^, which is introduced as 
part of the nonlinear gauge fixing condition. Similar to the maximal abelian gauge of quantum 
chromodynamics, the renormalization of the gauge parameter is singular. 
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1 Introduction. 



Gauge theories are the underlying quantum field theory describing the physics of the funda- 
mental properties of nature. For instance, the field theory describing the strong interaction, 
quantum chromodynamics (QCD), is based on a non-abelian colour group and exhibits asymp- 
totic freedom, [Il[2], whereby the basic quarks and gluons effectively behave as free particles at 
very high energies. Indeed this property allows one to apply perturbation theory to describe 
strong physics phenomena such as deep inelastic scattering to very high precision. One fun- 
damental reason for the basic feature of asymptotic freedom, [U [2], is the self-interaction of 
the gluon which is a natural consequence of the generalization of the abelian gauge theory of 
quantum electrodynamics (QED) to the non-abelian gauge group. The consequent nonlinearities 
introduced in this extension provide a richer structure not present in the dynamics of electrons 
and photons. In order to perform quantum calculations to determine the properties of a gauge 
theory, one has to first fix the gauge to ensure that only the physical degrees of freedom are taken 
into account. There is a large range of choices as to how to eliminate the unphysical degrees 
of freedom. In general such choices can be classified roughly under various headings, such as 
covariant or non-covariant, linear or nonlinear and physical or unphysical. (For a comprehensive 
review see, for example, [3J.) Ordinarily when one wishes to perform high precision and there- 
fore high loop calculations, one chooses linear covariant gauges (which are not physical) such as 
the Landau or Feynman gauge. This is primarily because their covariant, though unphysical, 
nature does not overcomplicate the resultant (massless) Feynman diagrams. By contrast other 
gauges, whilst being motivated by physical considerations, such as the Coulomb gauge, are not 
necessarily renormalizable. If they can be proved to be renormalizable at a formal level, then 
it is not always clear whether loop integrals beyond one loop can be computed, [3]. However, 
irrespective of how one fixes the gauge, one fundamental feature is always present and that is 
that physical predictions must always be independent of the choice of gauge. As the bulk of 
multiloop calculations are performed in linear covariant gauges and mass independent renor- 
malization schemes, one check in this instance is that the covariant gauge fixing parameter, a, 
must be absent in the determination of physical quantities or objects which are clearly gauge 
independent or gauge invariant, t4j. Indeed this property provides a powerful checking tool in 
intricate multiloop calculations. 

Although linear covariant gauges have been examined in detail, there has been recent interest 
in nonlinear covariant gauges due to their potential connection with the infrared dynamics of 
non-abelian gauge theories, [5l [6]. For instance, the Curci- Ferrari gauge, [7j, and maximal 
abelian gauge (MAG), [Hl[9l[T0], have been studied where the aim was to examine the abelian 
dominance hypothesis, [El [lU [l2l [13]. Briefly the gluons in the centre of the colour group are 
believed to dominate the infrared sector due to the off-diagonal gluons acquiring a dynamical 
mass greater than that of the centre gluons. The latter are then central to abelian monopole 
condensation that is believed to drive the confinement mechanism, [HJ \TT\ [T2| [T3] . However, 
from a practical point of view to perform any calculations in such nonlinear gauges one needs 
to renormalize the gauge theory to as high a loop order as is possible. This has been achieved 
for both the Curci- Ferrari gauge and MAG in the MS scheme at three loops, [14[ I15j. From 
the point of view of trying to understand basic features of covariant nonlinearly gauge fixed 
gauge theories both these gauges have common properties. For instance, unlike linear covariant 
gauges they have quartic ghost interactions and the corresponding gauge parameters get non- 
trivial renormalization. Further in the case of the MAG this gauge parameter renormalization 
is singular. Though gauge independent quantities, such as the /3-function, correctly emerge as 
gauge parameter independent. Whilst these gauges are essentially related by construction to 
the non-abelian aspect of the gauge theory itself, one natural question to ask is, is this a feature 
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in other nonlinear gauges when one has an abeUan structure. 

Clearly QED is invariably treated in a linear covariant gauge. However, with the explosion 
of interest in gauge theories in the early 1970's QED was gauge fixed in a nonlinear gauge known 
as the 't Hooft-Veltman gauge, [16], and shown to be renormalizable. The primary interest in 
this gauge fixing was that the nonlinearity naturally introduced an abelian gauge theory which 
mimicked QCD. This was due not only to the presence of interacting Faddeev-Popov ghosts but 
also triple and quartic photon self-interactions. As the latter do not appear in linear covariant 
gauges in QED, the 't Hooft-Veltman gauge clearly could be used as a laboratory to study simple 
issues related to gauge field self-interaction. Indeed in [T71 [18] a one loop calculation showed 
that the corresponding covariant gauge fixing parameter was renormalized. Therefore, in light 
of these observations it is the purpose of this article to record the full three loop renormalization 
of QED in the 't Hooft-Veltman gauge. Indeed as far as we are aware this will represent the first 
detailed multiloop study of the renormalization of QED in the 't Hooft-Veltman gauge. We will 
construct all the renormalization group functions in the MS scheme including the /3-function of 
the electron-photon coupling constant which will agree with the already established results of 

The paper is organised as follows. The relevant properties of QED gauge fixed in the 't Hooft- 
Veltman gauge are discussed in section two with the details of the full renormalization given in 
section three. Concluding comments are provided in section four. 



2 Background. 

First, we introduce the 't Hooft-Veltman gauge in QED, [16], and the notation and conventions 
we will use. The key ingredient is the gauge fixing functional, JF[^^], which slots into the 
conventional path integral formalism for constructing a quantized gauge theory. Here is the 
photon field. We take 

J^[A^] = d'^A^ + ^^Af'A^ (2.1) 

which is clearly nonlinear where for the moment ^ is a parameter and a is the gauge fixing 
parameter. Clearly when ^ = one recovers the usual linear gauge fixing functional whence a 
becomes equivalent to the gauge parameter of those gauges. In other studies of the 't Hooft- 
Veltman gauge, however, was invariably fixed to certain numerical values such as 1 or 2. We 
leave it as a free parameter here and given that eventually it will appear with the triple and 
quartic photon self-interactions we will regard it as a coupling constant which will run. It is 
not to be confused with the usual gauge coupling constant, e, which is present in the covariant 
derivative when electrons are present. Therefore we are in effect working with a two coupling 
theory. Though in the absence of electrons, whilst photon self-interactions are present the field 
theory is effectively a free theory of photons since the physics cannot be altered by the gauge 
fixing. This feature ought to emerge in the computations. Hence the full Lagrangian for Nf 
massless electrons in the 't Hooft-Veltman gauge is, |16j, 

L = - ^F^^'F^, + cd^'d^c + ^cA'^df.c + b (^df^A^ + ^C^^^/.) + ^«^' + ^V^W (2-2) 

where tp is the electron field, c and c are the Faddeev-Popov ghosts emerging from the path 
integral formalism and b is the Nakanishi-Lautrup auxiliary field which arises in the off-shell 
BRST formalism. Eliminating it by its equation of motion produces the Lagrangian in the form 
we will treat it 

L = - Jf'^^F^. + cd^'d^c + (cA^'d^c - ^ (o^'A^ + ^^A^'aX + iV^^ • (2-3) 
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Ordinarily in a linear covariant gauge in QED one drops the Faddeev-Popov ghosts from p.3p 
when = since they do not couple to photons or electrons. For ^ 7^ this is not possible 
and they are not only present but play a key role in the full renormalization of the theory. The 
covariant derivative, is defined by 

= + ieA^ . (2.4) 

Unlike the Curci-Ferrari gauge and the MAG in QCD, there is no quartic ghost self-interaction 
which is due to the absence of a colour index on the ghost fields meaning that c(x)c(x) = due 
to their anticommuting property. The Feynman rules for (|2.3p are straightforward to derive but 
the interested reader can view them in [23]. By construction (12. 2p is invariant under the BRST 
symmetry, [ig, 

6A^ = — dfj,c , 6c = , 6c = b , 6b = , 6'il> = iecip (2.5) 
which is clearly nilpotent. 



3 Renormalization. 



We turn now to the details of our three loop MS renormalization. If we regard the fields 
and parameters of (|2.3p as bare then to renormalize (12. 3p we introduce renormalized fields and 
variables by 
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where the subscript o denotes a bare quantity. We have also chosen to follow the same convention 
as [Hj for the definition of the renormalization of the gauge parameter a. With this convention 
it is the combination ZaZ~^ which is unity in the linear covariant gauge after renormalization. 
As we will use dimensional regularization in d = 4 — 2e dimensions with e as the regularizing 
parameter, the renormalization scale /U has been introduced to ensure that both renormalized 
couplings e and ^ remain dimensionless in d-dimensions. In principle all the renormalization 
constants will be functions of both (renormalized) couplings. However, it transpires that both 
the /3-function of e and the photon anomalous dimension are independent of ^. Therefore, the 
information determined from the MS renormalization will be encoded in the renormalization 
group functions. The /3- functions are given by 
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The anomalous dimensions are defined in the usual way by, (see, for example, [24]), 
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from which it is straightforward to deduce 
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in terms of the renormalization constants where we have set a = /{IQtt'^) and z = .^^/(IGvr^). 
In these definitions, (j3.4p . we have not assumed that 7^(0,2;) = 1 which is the case in a hnear 
covariant gauge. Clearly the wave function renormalizations will also depend on a. Moreover 
in these definitions where our explicit renormalization constants are clearly independent of one 
of the coupling constants or gauge parameter, we have included this property in the derivation 
of (j3.2p and (j3.4p . For instance, Zg turns out to independent of both a and a. 
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Two loop 


Three loop 


Total 
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Total 


9 


96 
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Table 1. Number of Feynman diagrams for the renormalization of each Green's function. 



The full three loop renormalization is performed for the massless case using the Mincer 
algorithm, |24l I25j . written in the symbolic manipulation language Form, |26]. The Feynman 
diagrams are generated automatically with the Qgraf package, [27], before being converted into 
Form input notation by converter routines. The number of Feynman diagrams we compute at 
each loop order for the set of Green's functions we need to consider to render (j2.3p finite are given 
in Table 1. The wave function renormalizations are deduced from the photon, Faddeev- Popov 
and electron 2-point functions whilst the 3-point functions determine the coupling constant 
renormalizations. For the latter given that the coupling constant renormalizations are gauge 
independent we have a strong check on the wave function calculations. Also for these, to apply 
the Mincer algorithm an external momentum has to be nullified. This is because Mincer 
computes massless 2-point functions up to the finite part at three loops. The extraction of 
each of the renormalization constants is found by applying the approach of [28] for automatic 
Feynman diagram calculations. The 2 or 3-point functions are computed as a function of the 
bare parameters. Then these are replaced by the renormalized variables from (j3.ip and the 
undetermined renormalization constant for that 2 or 3-point function chosen so as to absorb the 
infinities which remain. The latter appear as poles in e and are absorbed into the renormalization 
constants with the usual MS definition. Prior to presenting the results of our labours we note that 
one main check is that the double pole in e at two loops and the double and triple poles at three 
loops for any renormalization constant are predetermined by the structure of the renormalization 
group equation. In the expressions we present for the anomalous dimensions and /3-functions all 
the renormalization constants passed this test. 

Hence, the complete set of three loop MS renormalization group functions are 
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where the formal order symbol 0(z"a'~") means all appropriate possible combinations of the 
coupling constants a and z at the Ith loop. In both /3-functions the d-dimensional dependence has 
been retained as an indication of our conventions in deriving the renormalization group functions 
as well as for the reader interested in constructing the original renormalization constants from 
the differential equations of (|3.2p and (j3.4p . 

Aside from the internal checks based on consistency with the renormalization group equation 
there are additional checks on these results. First, both /3-functions correctly emerge as a 
independent and f3a{a) is in total agreement with the original linear covariant gauge result of 
[191 l20l [2H I22j . Also the photon anomalous dimension is proportional to [3a{a) as required by 
the Ward identity. Further, in the limit z ^ the anomalous dimensions agree with those of 
the linear covariant gauge fixing, [28]. Finally, one must recover a theory of free photons when 
the electron interaction is switched off via a ^ 0. Clearly, /^^(O, z) = which corresponds to 
a free field theory when ^ 7^ even though neither 7^(0,2:), 7c(0, z) nor 7^(0, z) are zero. An 
additional comment on our results is that 7^(0, 0) has clearly only a dependence at one loop. 
This was originally observed in [29], where it was claimed that the only a dependence of 7^(0, 0) 
was at one loop. Clearly in a nonlinear covariant gauge there is a dependence beyond one loop 
which is not unexpected. Finally, in relation to the renormalization of the gauge parameter in 
other conventions, we record the sum of 7a(«) and 7a(a, z) is 
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which is clearly non-zero for z ^ 0. Moreover, like the MAG, (see, for instance, [IS]), the 
corresponding anomalous dimension is also singular in the a — > limit, though similarly the 
remaining renormalization group functions, including the /3-functions, are finite in this limit. 
For 7A(a), 7c(a, 2) and ^^{a,z) this is primarily because in (j3.3|) and (j3.4p . the term involving 
7Q,(a, z) is multiplied by a and Za, Zc and Z^ themselves are non-singular at a = 0. 

Finally, it might be tempting to try and remove the a = singularity in 70(0, z) by a suit- 
able coupling constant redefinition. Whilst this would produce renormalization group functions 
analytic in a, one must be careful in ensuring that the original theory is retained. For instance, 
to remove the 1/a terms in 70,(0, z) the least one must do is to redefine z by a factor proportional 
to a. In the simplest case, setting z = az one would formally have non-singular renormaliza- 
tion group functions. However, in this instance returning to (j2.3p in the absence of electrons 
the Landau gauge Lagrangian would then describe self-interacting photons with non-interacting 
ghosts. This is not consistent with the notion that without electrons the photon is a free field. 
Therefore given this, avoiding what might be perceived to be a problem in a renormalization 
group function, which has no physical interpretation, in order to render it analytic, has a signif- 
icant affect on the nature of the original theory. In other words whilst it might seem unnatural 
to have a theory with couplings which are singular as a — > leading to a singular anomalous 
dimension, the nature of the theory remains consistent. 



4 Discussion. 



We have completed the full three loop renormalization of QED in the nonlinear 't Hooft-Veltman 
gauge. This extends the one loop calculations of [17], HI]. Whilst the authors of [17], [18] were the 
first to observe that the longitudinal part of the photon is renormalized unlike in a linear gauge, 
we have carried out a slightly more general analysis by allowing for a covariant gauge parameter 
a and the inclusion of an additional coupling ^ in order to track the loop calculation in a similar 
way to the usual coupling constant a. The coupling ^ was not initially fixed to a specific value. 
Consequently the singular renormalization of the gauge parameter emerges. Whilst this is not a 
new feature of a nonlinear gauge fixing, since the MAG of QCD has the same property, it does 
not disrupt either the renormalizabilty of the theory, [16] , or the evaluation of gauge independent 
quantities such as the /3-functions. This is primarily because although the gauge in one sense 
is only defined in the a — > limit, (j2.3p . the anomalous dimension of a, ()3.6p . has no physical 
meaning or interpretation. 
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